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, Abstract 

OO 

OO ■ We use detailed balance for a hadron composed of quark and gluon 

, Fock states to obtain parton distributions in the proton and pion on 

I the basis of a simple statistical model. 
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g^' 1 Introduction: parton distributions in the 



proton 



^ ■ There has been considerable interest in the flavor dependence of the proton's 

. quark and antiquark distributions. The first measurement of the {d — u) 

asymmetry was made by the NMC group The integral of this distribution 
showed a violation of the Gottfried sum rule. Later, Drell-Yan |2j||3j and deep 
inelastic scattering experiments determined the Bjorken-x dependence of 
the asymmetry. The meson cloud model and the Sullivan process were used 
to explain the momentum fraction distribution of {d — u); many other models 
have been proposed [S]. 

Most recently, Zhang and collaborators [^13111 have used a simple statis- 
tical model to calculate the (d — u) distribution in the proton. They consider 
the proton to be an ensemble of quark-gluon Fock states, and use detailed 
balance 13 o^' "the principle of balance" jH] to determine the distribution 
functions for all partons of the proton. Despite its simplicity, the model does 
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reasonably well in predicting the distributions of partons, as well as for that 
of {d — u). The excess of d over u comes about from the 2:1 ratio of u:d, 
which provides an excess of u quarks for the annihilation of m's. 

Zhang, Zou, and Yang (ZZY) J_j write a general Fock state expansion for 
the proton as 

\P >= J2 Cijk\{uud}{ijk} >, (1) 

with i the number of uu pairs, j the number of dd pairs and k the number 
of gluons. The states are normalized such that the sum of the probabilities 
Pijk = |cjjfcP of finding a proton in the state \{uud}{ijk} > summed over all 
i, j, and k is unity 0, 

E P^J^ = 1 • (2) 

In ZZY's statistical model, detailed balance between any two Fock states 
requires that 



PijkN{\{uud}{ijk} >^ \{uud}{i'j'k'} >) = 
Pi>j>k'Ni\{uud}{i'j'k'} >^ \{uud}{ijk}) >, (3) 

in which N{A ^ B)) is the transfer rate of state A into state B. Transfer 
rates between states are assumed to be proportional to the number of partons 
that can split or recombine. Taking into account two processes, q ^ q g and 
g ^ qq, ZZY find that 

EMI = I (A) 

pooo tKt + 2)\j\ij + l)\k\- 

This equation, together with the normalization condition determines all 
the Pijk- It is clear from this equation that uu states, labelled by i, are 
suppressed relative to dd states, labelled by j. Summing over all states, 
Zhang et al. jH] find {d — u) ~ 0.124, remarkably close to the experimental 
value of 0.118 ±0.012 g]. 

ZZY determined parton distribution functions for the proton by using a 
Monte Carlo simulation of the distribution of momenta among the n par- 
tons in each Fock state. The phase space volume for n free partons is 
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determined by 

c^^v-gponpl^. (5) 

with P and pi the 4-momenta of the proton and the z-th parton, respectively. 
The masses of the partons are neglected so that = and 

n P rl P rIO 

c=^^(i^-E/'.)n^^^^^. (6) 

ZZY argue that this free parton phase space distribution should be multi- 
plied by n -ET^ because partons with smaller momenta spend more time at 
the center of the proton where they are almost free; these partons are thus 
weighted with a higher probability. Then dfn, the distribution for confined 
partons, is: 

We have found that the effect of the weighting factor is quite small, except 
for the very lowest and highest parton momenta. 

From the Monte Carlo distribution of parton momenta Pi, the parton 
distributions can be found in terms of the light cone variable Bjorken x, 

Ei Pzi 

Xi = — — , (8) 

in which M is the proton mass. We used RAMBO [H] for our Monte Carlo 
event generator. Then for an ?7,-parton state, for which n = 3 + 2{i + j) + k, 
the ^-distributions for u and d are 



Uijk{x) = fn{x)i , dijk{x) = fn{x)j , (9) 

for u and d are 

Uijk{x) = fn{x){2 + i) , dijk{x) = fn{x){l + j) , (10) 
and for the gluons is 

gijk{x) = fn{x)k . (11) 
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Thus, we find, in accord with ZZY, 

Uix) = J2 PijkUijk{x) , (12) 

and corresponding equations for d{x), u{x), d{x) and g{x), normalized so 
that 

/ x[u{x) + d{x) + u{x) + d{x) + g{x)\dx = 1 . (13) 
The average number of partons in the proton, ra, is given by 

n = [u{x) + d{x) + u(x) + d{x) + g{x)]dx ^ 5.6 . (14) 
Jo 

This sets the scale /io ~ -E = M/n ^0.17 GeV for Ql = at which the 
distributions are calculated. The plot in Fig. 1, which reproduces ZZY's 
results, shows that the experimentally deduced (J — u) is fit qualitatively in 
this model , but is low at small x and high at large x. The discrepancy with 
experiment shows up more starkly in Fig. 2, our plot of d{x)/u{x). We find 
that these results are changed very little if the phase space weighting factors 
O-ET™ varied from m = 1 to m = 3. Nevertheless, we believe that it is 
remarkable that such a simple model does so well. 

There are, of course, other explanations of the excess of d over u, which 
have been considered previously, particularly the pion cloud [TUl HH US [131 
ITl] and uj mesons The leading term in the Fock state expansion of the 
pion cloud model is a "bare" proton which consists of valence quarks plus 
qq pairs due to gluon splitting. It can be argued that this gluon splitting is 
what is being considered by ZZY and thus the perturbative sea in the "bare" 
proton should not be symmetric, i.e. d ^ u. The pion cloud, represented by 
the higher-order terms in the expansion, would then be an additional effect. 
If this is the case, the pions must have a larger role at high x where the 
fall-off of the ratio d/u is not reproduced by the statistical model. We will 
not pursue this argument further here. 

2 Parton distributions in the pion 

If the statistical model has some validity, then it should not only work for the 
proton, but also for the pion. This distribution function is, in fact, required 
in the pion cloud model. We have therefore investigated the valence and sea 
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Figure 1: Comparison of statistical model calculation with E866 experimental 
results 12] for J — M. 



5 



1.5 - 



d/u 



0.5 - 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 



Figure 2: Comparison of statistical model calculation with E866 experimental 
results 12] for d/u. 
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quark distributions for the tt"*" in the statistical model. The formulas are 
similar to those for the proton, but because there is only one valence quark 
of each flavor in the tt"*", the sea is flavor symmetric. 
We write the Fock state expansion for the pion as 

k+ >= Cijk\{ud}{ijk} > . (15) 

The analysis of section 1 is unchanged, except that n = 2 + 2{i + j) + k, and 
the ratio of probabilities for different Fock states is now 

so that the Tr"*" sea is symmetric, i.e. u{x) = d{x). We find n = 4.5 in the 
pion. 

The quark distributions for a specific n-parton state are 

dijk{x) = fn{x){l +j) , Uijk{x) = /n(x)(l + i) , (17) 



Uijk{x) = fn{x)i , dijk{x) = fn{x)j , (18) 

and for the gluons is 

gijk{x) = fn{x)k . (19) 

The parton distributions are found by summing these distribution functions 
over all values of {ijk} 

= PijkUijk{x) = d{x) , (20) 
d{x) = Y Pijk dijk{x) = u{x) , (21) 

and 

gix) = Y Pijk 9ijkix) . (22) 

The valence quark distribution function is 

v{x) = u{x) — u{x) = d{x) — d{x) (23) 
Our results are shown in Fig. 3. The valence quark distributions are too high 
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Figure 3: Our results for parton density distributions xq{x) and xg{x) for 
the pion. Solid curve: valence quark distribution; long-dashed curve: sea 
quark distribution; short-dashed curve: gluon distribution. 
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Figure 4: Our calculation of the valence quark distribution x v{x) in the pion, 
compared to the experimental results of Conway et al. |W . The dashed curve 
shows our results without any evolution, as in Fig. 3, which correspond to 
a scale of Ql = 1.96 GeV^. The solid curve shows our results evolved to 
Q"^ = 16 GeV^ of the E615 experiment. 
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for large x because of the dominant contribution of the n = 2, {ijk} = {000} 
state, the leading term in the Fock expansion, for which f2{x) is a uniform 
distribution in x. The sea quark distribution is flavor symmetric, as noted 
above. In order to compare our valence quark distributions to experiment, 
we carried out an evolution in Q^. We determined the starting scale of our 
distributions by requiring that the first and second moments of our valence 
quark distribution at = 4 GeV^ be equal to those found by Sutton et al. 
[T7j. This gave us a starting scale of = 1.96 GeV^. We used Miyama and 
Kumano's code BFl [18j for the DGLAP evolution. We compare our pion 
valence quark distribution with that obtained by E615 from pion scattering 
on tungsten [in] in Fig. 4. The dashed curve shows our results without 
any evolution, as in Fig. 3. The solid curve shows our results evolved to 
= 16 GeV^ of the E615 experiment. The agreement between theory and 
experiment is good. Other theoretical calculations of pion parton distribu- 
tion functions have used constituent quark models, the Nambu-Jona-Lasinio 
model, instantons, or the Dyson Schwinger equations. Moments of the dis- 
tributions can be calculated in lattice QCD, from which particular forms of 
the distributions can be reconstructed. For references see the recent papers 
of Hecht, Roberts and Schmidt [20] and Detmold, Melnitchouk and Thomas 
pij . In Fig. 5 we compare our valence quark distribution to the Dyson- 
Schwinger calculation of Hecht et al. and to experiment. Both distributions 
were evolved to = 16 GeV^. We find it remarkable that our simple sta- 
tistical model agrees with experiment as well as the covariant, QCD-based 
model. 



3 Conclusions 

The calculation of parton distribution functions is an important goal of non- 
perturbative QCD. We have used the statistical model of Zhang et al., de- 
veloped for the calculation of parton distribution functions in the proton, 
to calculate the parton distribution functions of the pion. We find that this 
simple model, with no free parameters, is in good agreement with experiment 
and other calculations. 
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Figure 5: Our results (solid curve) for the valence quark distribution xv{x) 
in the pion, compared to the calculation of Hecht, Roberts and Smith 2(]] 
(dashed curve) and the experimental results of Conway et al. [16 . Both 
calculations were evolved to = 16 GeV^ of the E615 experiment. 



11 



4 Acknowledgement 



We thank Toby Burnett and Steve Ellis for discussions on the use of RAMBO, 
and Michael Clement for computational work with RAMBO and calculation 
of distribution functions. This work was supported in part by the U.S. De- 
partment of Energy and the Research in Undergraduate Institutions program 
of the U.S. National Science Foundation, Grant No.0070942. 

References 

[1] R Amaudruz et al., NMC Collaboration, Rhys. Rev. Lett. 66 (1991) 
2712. 

[2] A. Baldit et al., Rhys. Lett. 332 (1994) 244. 

[3] E866 Collaboration, E.A. Hawker et al., Phys. Rev. Lett. 80 (1998) 
3715 ; J.C. Peng et al., Phys. Rev D58 (1999) 092004; R.S. Towell et 
al., Phys. Rev. D64 (2001) 052002. 

[4] K. Ackerstaff, et al, Phys. Rev. Lett. 81(1998) 5519. 

[5] J. Speth and A.W. Thomas, Advances in Nuclear Physics, Vol. 24, ed. 
J.W. Negele and E.W. Vogt (Plenum Press, NY) 1998; S. Kumano, 
Phys. Rep. 303 (1998) 183. 

[6] Y.-J. Zhang, B. Zhang and B.-Q. Ma, Phys. Lett. B 523 (2001) 260. 

[7] Y.-J. Zhang, B.-S. Zou and L.-M. Yang, Phys. Lett. B 528 (2002) 228. 

[8] Y.-J. Zhang, W.-Z. Deng and B.-Q. Ma, Phys. Rev. D 65 (2002) 114005. 

[9] R. Kleiss, W. J. Stirhng and S. D. Ellis, Comp. Phys. Comm. 40 (1986) 
359. 

[10] S. Kumano, Phys. Rev. D43 (1991) 3067. 

[11] A. Signal, A.W. Schreiber and A.W. Thomas, Mod. Phys. Lett. A6 
(1991) 271. 

[12] W. Koepf, L.L. Frankfurt and M. Strikman, Phys. Rev. D53 (1996) 
2586. 



12 



[13] W. Melnitchouk, J. Speth, and A.W. Thomas, Phys. Rev. D59 (1998) 
014033. 

[14] N.N. Nikolaev et al., Phys. Rev. D 60 (1999) 014004. 

[15] M. Alberg, E.M. Henley and G.A. Miller, Phys. Lett. B471 (2000)396. 

[16] J.S. Conway et al., Phys. Rev. D 39 (1989) 92. 

[17] P.J. Sutton et al., Phys. Rev. D 45 (1992) 2349. 

[18] M. Miyama and S. Kumano, Comp. Phys. Comm. 94 (1996) 185. 

[19] V.N. Gribov and L.N. Lipatov, Sov. J. Nucl. Phys. 15 (1972) 438 and 
675; G. Altarclh and G. Parisi, Nucl. Phys. B 126 (1977) 298; Y. L. 
Dokshitzer, Sov. Phys. JETP 46 (1977) 641. 

[20] M.B. Hecht, CD. Roberts and S. M. Schmidt, Phys. Rev. C 63 (2001) 
025213. 

[21] W. Detmold, W. Melnitchouk and A. W. Thomas, Phys. Rev. D 68 
(2003) 034025. 



13 



